arXiv: 1503.00567v2 [cond-mat.mtrl-sci] 14 May 2015 


A many—body perturbation theory approach to the electron—phonon interaction with 

density—functional theory as a starting point 


Andrea Marini, 1,2 S. Ponce, 3,2 and X. Gonze 3,2 

1 Istituto di Struttura della Materia of the National Research Council, 

Via Salaria Km 29.3, 1-00016 Monterotondo Stazione, Italy 
2 European Theoretical Spectroscopy Facilities (ETSF) 

3 Universite catholique de Louvain, Institute of Condensed Matter and Nanosciences, 

NAPS Chemin des etoiles 8, bte L07.03.01, B-13^8 Louvain-la-neuve, Belgium 

(PDated: May 15, 2015) 

The electron-phonon interaction plays a crucial role in many fields of physics and chemistry. 
Nevertheless, its actual calculation by means of modern many-body perturbation theory is weak¬ 
ened by the use of model Hamiltonians that are based on parameters difficult to extract from the 
experiments. Such shortcoming can be bypassed by using density-functional theory to evaluate 
the electron-phonon scattering amplitudes, phonon frequencies and electronic bare energies. In this 
work, we discuss how a consistent many-body diagrammatic expansion can be constructed on top of 
density-functional theory. In that context, the role played by screening and self-consistency when 
all the components of the electron-nucleus and nucleus-nucleus interactions are taken into account 
is paramount. A way to avoid overscreening is notably presented. Finally, we derive cancellations 
rules as well as internal consistency constraints in order to draw a clear, sound and practical scheme 
to merge many-body perturbation and density-functional theory. 

PACS numbers: 71.38.-k,71.15.Mb, 71.15.-m 


In physics and chemistry the interaction between elec¬ 
tronic and vibrational degrees of freedom is at the ori¬ 
gin of a multitude of phenomena. Focusing on solid 
state physics, this coupling usually determines the elec¬ 
trical and thermal conductivity of metals as well as car¬ 
rier lifetime in doped semiconductors^. It also induces 
the transition to a superconducting phase in many solids 
and nanostructures^. The electron-nucleus coupling also 
plays a role in the renormalization of electronic bands^, 
carriers mobility in organic devices^ and dissociation at 
the donor/acceptor interface in organic photovoltaic^. 
This coupling can naturally interact with other couplings 
like the magnetic field leading, for example, to the spin- 
Seebeck effect^. 

This effect is nowadays the subject of an intense re¬ 
search activity for its crucial role in new emerging fields 
of experimental and theoretical physics. The electron 
nucleus coupling plays a crucial role in the relaxation 
and dissipation of photo-excited carriers in pump and 
probe experiments^. Similarly, modern Angle-Resolved 
Photoemission Spectra (ARPES) experiments have re¬ 
cently disclosed the complex and temperature depen¬ 
dent structures appearing in the spectral functions of 
several oxides^. These structures quite remarkably re¬ 
semble similar structures predicted to exist in conjugated 
polymer s 10 : 11 pointing to a strong effect whose physical 
origin is still not completely clear. 

From the theoretical point of view, the most up- 
to-date scheme to calculate and predict the ground- 
and excited-state properties of a wide range of mate¬ 
rials is based on the merging of Density-Functional- 
Theory (DFTU2 with Many-Body Perturbation The¬ 
ory (MBPT)±2. 

DFT is a broadly used ab-initio ground-state theory, 


that allows to calculate exactly electronic density and 
total energy without adjustable parameter. The merg¬ 
ing of DFT with perturbation theory gives the so-called 
Density-Functional Perturbation TheoryiA— (DFPT). 
The DFPT is a powerful computational tool for the direct 
treatment of phonons. 

However, the DFT computation of excited electronic 
states properties like the bandgap energies is a known 
problematic topioi^. As a result, MBPT is nowadays the 
preferred alternative to DFT for that purpose. It is based 
on the accurate treatment of correlation effects by means 
of the Green’s function formalism. MBPT is formally 
correct and leads to a close agreement with experiment— i 
but is extremely computationally demanding. A natu¬ 
ral way to solve this issue is to merge the quick DFT 
calculation with the accurate MBPT one. The latter 
method is often referred to as ab-initio Many-Body Per¬ 
turbation Theory^ (ai-MBPT). In this method, DFT 
provides a suitable single-particle basis for the MBPT 
scheme. This methods has been applied successfully to 
correct the well-known band-gap underestimation prob¬ 
lem of DFT-i&ii. 

Although the ai-MBPT aims at calculating the excited 
state properties with an unprecedent precision, it is com¬ 
monly applied by neglecting the effect of lattice vibra¬ 
tions. Even today, most of the ai-MBPT results are com¬ 
pared with finite-temperature experimental data—. Such 
comparison is not even well motivated at zero tempera¬ 
ture as the lattice vibrations induce a zero-point motion 
effect that can be sizeable— i 11 ! 21 ! 22 , e.g. on the order of 
0.4-0.6 eV for the direct and indirect band gaps of dia- 
mon d 23 ! 24 . This represent a clear motivation to develop a 
coherent ab-initio theory in which the electron-phonon 
interaction is rigorously included on top of ai-MBPT. 
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The need for such theory is exemplified by the very 
fragmented historical development of the ab-initio ap¬ 
proach to the temperature dependence of the electronic 
structure due to the electron-phonon interaction. From 
the fifties to the late eighties, a coherent ab-initio frame¬ 
work was still not devised and the electron-phonon in¬ 
teraction was initially investigated and computed in a 
semi-empirical context by Fa n 25 i 26 . His theory had no 
adjustable parameters and was based on the first-order 
perturbed Hamiltonian. During the same period, An- 
toncil—, followed by others^Sr— , developed empirical 
Debye-Waller (DW) corrections to the nuclear potential. 
Only in 1976, Allen and Heine^i rigorously unified the 
Fan and DW corrections in a common framework. Their 
approach, combined with the use of a semi-empirical 
model, allows for a re-writing of the problem in terms of 
first-order derivatives of an effective potential only. Cal¬ 
culations of the electron-phonon renormalization effects 
were then led by Cardona and coworkers2i22ri2i, including 
Allen. The resulting approach is now called the Allcn- 
Heine-Cardona (AHC) theory. 

In 1989, the first ab-initio calculation of the temper¬ 
ature dependence of the gap was attempted, by King- 
Smith et al — , based on DFPT. Starting from there, sev¬ 
eral first-principle calculations have been done, relying 
mainly on three types of formalisms: (i) time averaging 
of bandgap obtained using first principles molecular dy¬ 
namics simulations^ - —; (ii) frozen phonons (FP) calcu¬ 
lation s 24 - 41- — and (iii) the AHC approach implemented 
in a full ab-initio framework by using DFT and DFPT 
as a reference system 2 - 2 -24 : 45 ! 46 . 

All these approaches are based on an adiabatic and 
static treatment of the electron-phonon interaction. This 
limitation was overcome by using the dynamical version 
of MBPT by Marini et a / 10 ' 11 ' 47 who focused on retarda¬ 
tion effects. 

Since then, there have been an increasing number of 
studies in which the electron-phonon interaction is fully 
included in the computation of the electronic structure, 
well beyond DFT. Still, several basic questions remain. 
In particular, the use of an electron-phonon interaction 
whose strength is computed from DFT, in a formalism 
that goes beyond DFT, e.g. the ai-MBPT approach, 
leads to several ambiguities due to the simultaneous in¬ 
clusion of different levels of correlation at the MBPT and 
DFT/DFPT level. 

Indeed, the DFPT electron-phonon interaction is nat¬ 
urally screened as it is computed from the derivative of 
the self-consistent Kohn-Sham potential with respect to 
atomic displacements. This screening is taken as it is 
in the MBPT part of the a*-MBPT scheme, although 
it is well known that the diagrammatic technique also 
predicts the screening of the electron-phonon interaction 
consistently with the kind of correlation included in the 
self-energ y 48 ! 49 . It has in fact been shown that the size 
of the zero-point motion renormalization is significantly 
larger in the MBPT than in the DFPT approaches— 

Another important issue of the az-MBPT approach is 


the lack of a diagrammatic interpretation of the screen¬ 
ing of the Debye-Waller term. This screening arises quite 
naturally in the DFPT— and AHC approaches— It is 
easy to show that it comes from the DFT self-consistent 
screened ionic potential. Instead, in the pure MBPT 
treatment of electrons and nuclei, this diagram is un¬ 
screened. The Debye-Waller diagram is however usually 
taken as screened without justification in most practical 
application because of the DFPT basis. 

In addition, a non-rigid nuclei correction to the Debye- 
Waller contribution^ is predicted to exist within the 
DFPT approach. However, this term is notably absent 
from the standard derivation of the electron-phonon the¬ 
ory based on the MBPT. 

The last issue is even more fundamental. Most of the 
electron-phonon interaction treatments that appears in 
textbooks, see e.g. Refs. [48l - [50| . are based on the study 
of the homogeneous electronic gas (jellium). At variance 
with any realistic material, the jellium model is based 
on a drastic approximation: the ions are replaced by a 
jelly of positive charge, in contrast with realistic mate¬ 
rials where the nuclei and their mutual interaction must 
be taken explicitly into account. This is correctly done in 
DFT and DFPT but not in the MBPT approach derived 
from the jellium model. 

This paper aims at answering all these questions by de¬ 
vising a coherent, formal and accurate approach to merge 
the MBPT scheme with DFT and DFPT. We present a 
consistent electron-phonon interaction theory based on 
the MBPT formalism, insisting specifically on the con¬ 
nection between the MBPT and AHC approaches. This 
work is inspired by the seminal works of Allen^i and van 
Leeuwen^, going further by including the full description 
of the atomic potential into account. 

The merging procedure will lead to the natural defini¬ 
tion of a series of practical rules and advices about how 
to perform electron-phonon calculations on top of DFT 
without double counting problems. These series of rules 
are well justified within the ai-MBPT scheme that, in 
its practical form used in material science calculations, 
can be seen as a collection of prescriptions only partially 
based on a solid theoretical ground and rather inspired 
by the succesfull comparison with the experiments of sev¬ 
eral different materials. This a posteriori validation rep¬ 
resents and important part of the ai-MBPT approach. 

The structure of the paper is as follows. Section Q] 
presents the total Hamiltonian and introduces the no¬ 
tation. In section El we draw a parallel between the 
electron-electron and the electron-phonon self-energies 
to show what is the source of the problems that arise in 
the merging of MBPT with DFT and DFPT. Section Ell 
properly defines the reference Hamiltonian to be used as 
a zero-th order in the many-body expansion. In Sec¬ 
tion UVj the different interaction terms are described, in¬ 
cluding the contributions from the nuclei-nuclei interac¬ 
tion. In Section [V] we perform the formal diagrammatic 
summations at different level of approximations: Hartree, 
Hartree-Fock and GW. We use the different levels of cor- 
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relation of these self-energies to discuss the different role 
played by self-consistency diagrams and how the screen¬ 
ing of the interaction terms arises. At the same time we 
derive cancellation rules that highlight the crucial role 
played by the nuclei-nuclei interaction. 

Finally, Section fVTl reviews the DFPT approach to the 
electron-phonon coupling in order, in Section m to 
compare the different properties of the DFT and of the 
Many-Body approach. We provide, in a practical and 
schematic way, a series of formal properties of the Many- 
Body expansion performed on top of the DFT reference 
Hamiltonian. We discuss, from a diagrammatic perspec¬ 
tive, the physical origin of the Debye-Waller terms be¬ 
yond the screened rigid-ion contribution (Section IVII CD 
and a practical approach to calculate iteratively the n— 
th order derivatives of the DFT self-consistent potential 
(Section TYIIDl and Appendix[B]). Atomic (Hartree) units 
are used throughout the article. 

I. THE TOTAL HAMILTONIAN 


operator of the electron i and v (r — r') = |r — r'| 1 is 
the bare Coulomb potential. Similarly, 

W n -n (R) = i Y! ls Va , ZsZ ° ,V - Riv) > (5) 

W e-e=^E ( 6 ) 

with . = V. , 

We now use the notation O (R), or equivalently O (R), 
to indicate a quantity or an operator that is evaluated 
with the nuclei frozen in their equilibrium crystallo¬ 
graphic positions (R). We expand the Hamiltonian as 
a Taylor series up to second order in the nuclear dis¬ 
placements, 

H (R) « H (R) + Y, d Rlsa H (R)AR lsa 

Isa 

+ 2 dRi sa R l , B ,pH(R-)^R isa ^R i ' s 'Pi (?) 

lsa,l f s'/3 


We start from the generic form of the total Hamilto¬ 
nian of the system, that we divide in its electronic H e , nu¬ 
clear H n (R) and electron-nucleus (e-n) W e _ n (R) con¬ 
tributions, 


H (R) = H e + H n (R) + W e _ n (R), (1) 

where R is a generic notation that represents a depen¬ 
dence on the positions of the nuclei. 

The electronic and nuclear parts are divided in a ki¬ 
netic T and interaction part W : 

H e = T e + W e . e , (2) 

H n (R) =T n + W n _ n (R). (3) 

Note that the nuclear kinetic energy depends on the nu¬ 
clear momenta, and not on the actual positions of the 
nuclei. In the above definitions, the operators are bare 
(un-dressed). The analysis of the dressing of W e - n that 
arises as a consequence of the electronic correlations is 
one of the key objective of this work. Indeed, in the 
Many-Body (MB) approach, this dressing appears in the 
perturbative expansion in the form of electron-hole pair 
excitations and therefore, cannot be introduced a priori 
in the definition of the Hamiltonian. 

The explicit expression for the bare (e-n) interaction 
term is 


where a and /? are Cartesian coordinates and 

ARi sa = ^Risa — Ri sa lj . (8) 

The equilibrium crystallographic positions R are de¬ 
fined, in the present context, as the positions minimizing 
the expectation energy of the Born-Oppenheimer Hamil¬ 
tonian (with fully correlated electrons), i.e. all the con¬ 
tributions to the total Hamiltonian, except the nuclear 
kinetic energy, 

H BO (R) = H e + W n _ n (R) + W e _ n (R). (9) 


Those positions are equivalently defined by the condition 
that the expectation of the Born-Oppenheimer force Fr 1b 
acting on the nucleus located at position Rj s is zero 


Fr !s = -d Rlg (Hro (R.)) 


0 V{l,s}. (10) 


The average in Eq. m is done on the exact elec¬ 
tronic ground state of the Born-Oppenheimer Hamilto¬ 
nian. Still, the present theory will go beyond the Born- 
Oppenheimer approximation by considering fluctuations 
around the equilibrium positions. 


II. THE PROBLEM 


VF e _ n (R) = 



where Rj s is the nuclear position operator for the nucleus 
s inside the cell l (the cell is located at position R;), Z s 
is the corresponding charge, r, is the electronic position 


The problem we aim at solving is how to treat the ef¬ 
fect of the two last terms in the right-hand side of Eq. 0 
and how to do it by merging the MB approach, well- 
established for the treatment of Hro (R), with a DFT 
description of the reference electronic and nuclear sys¬ 
tems. 

When H (R) « H BO (R) the Hamiltonian represents 
indeed a purely electronic problem, for which the MB 
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approach is well-established in the literatur e) 13 ’ 49 ! 50 ' 53 . 
It relies on the definition of an electronic self-energy 
E(r,r';w) that is a complex and non local function in 
frequency and space. E can be approximated by follow¬ 
ing different strategies available in the literature (like the 
well known GW approximation^). For periodic solids, 
once the self-energy is known, the calculation of the cor¬ 
rection to an energy level |?ik) can be obtained by solving 
the corresponding Dyson equation (n is a band index and 
k the corresponding wavevector). 

Usually, the MB methodology starts from an 
independent-particle (IP) electronic Hamiltonian that 
includes only the kinetic electronic operator and the 
electron-nucleus operator, 

Hip (R) =T e + W e - n (R) . (11) 

The analysis of the correlated electronic Hamiltonian, 

H corr (R) = Hip (R) + W e _ e , (12) 

is addressed through the diagrammatic expansion. A 
simple approximation to the solution of the Dyson Equa¬ 
tion that fully captures the role played by correlation 
effects is the on-the-mass-shell approximation where: 

e«k « eik + 0 (™ k l s (u r '; e i°k) l nk )o, (13) 

where |nk)o and are the n-th single-particle eigen¬ 
state and eigenenergy of the independent-particle (IP) 
Hamiltonian Hip (R). 

In the present context, where we must consider dif¬ 
ferent configurations of nuclei, and determine also the 
equilibrium geometry through Eq. m, the initial cor¬ 
relation present in the reference system for the diagram¬ 
matic expansion must be carefully analyzed. Adding the 
nucleus-nucleus energy to Hip (R) gives 

H 0 (R) = T e + W e _ n (R) + W n _ n (R), (14) 

namely, the Born-Oppenheimer Hamiltonian without 
electron-electron interaction operator W e - e . This ini¬ 
tial Hamiltonian Hq do not include electron-electron cor¬ 
relations, and can be used as the starting point of a 
MB approach to the electronic problem. However, using 
this Hamiltonian instead of the true Born-Oppenheimer 
Hamiltonian means that no electronic correlation en¬ 
ergy contribution appears in the total energy and in the 
definition of the equilibrium nuclear positions through 
Eq. (USD. This would lead to a completely irrealistic de¬ 
scription of the starting nuclear geometry and vibrational 
frequencies. Indeed, e.g. the latter could be imaginary, 
and this would lead to unusual technical problems with 
the canonical transformation from the displacement op¬ 
erator to the phonon creation and annihilation operators. 
Thus, the extension of the electronic-only MB approach 
to the case where nuclear displacements are considered 
cannot use such a starting point. 


As DFT provides a treatment of energy and forces 
that include the electron-electron interaction, it yields 
a better starting point than Eq. (fill) . DFT is an exact 
mean-field theory in the sense that all electronic corre¬ 
lation effects are embodied in a mean-field exchange- 
correlation (xc) potential V xc [p] (?) which replaces the 
full electron-electron interaction operator W e - e , and de¬ 
pends on the ground-state density p. The bracket [] in 
V xc denotes a functional dependence. 

By adding to V xc \p\ (?) the Hartree potential, 
Vp [p\ (?), we get the total DFT potential: 

Vhxc [p] (?) =Y1 Vh M + V *c [P\ (?i), (15) 

i 

with 

V H [p] (?) = J dr'v (? - r') p (r'). (16) 

DFT is exact in the sense that the corresponding Kohn- 
Sham (KS) Hamiltonian 

Hks (R) = T e + V Hxc [p] + W e _ n (R), (17) 

provides, when the nuclear positions are given, a set of 
electronic eigenvectors whose corresponding density is 
the exact ground state density of Hpo (R) (Hohenberg- 
Kohn theorem^). 

The Hohenberg-Kohn theorem also states that Vh xc [p\ 
and the ground-state energy are functional of the exact 
electronic density p. It follows that, once the correct 
exchange-correlation functional is used, DFT gives the 
exact equilibrium nuclear positions through Eq. CO- 

In practice, an exact expression for Vhxc is not known 
and several approximations for it have been proposed 
in the literature^. In any case, even the simple local- 
density approximation (LDA) 55 ’ 56 , provides quite reason¬ 
able structural properties. Thus DFT represents a con¬ 
crete and accurate reference Hamiltonian to be used as 
zero-th order for a diagrammatic expansion that will al¬ 
low vibrational degrees of freedom to be included. For¬ 
mally, at the equilibrium geometry, one decomposes the 
correlated Hamiltonian as 

H corr (R) = H ks (R) + We_ e - V Hxc [p] - (18) 

At this point, the perturbative expansion is performed 
in terms of W e - e — Vh xc instead of W e - e . This is the 
theoretical basis of the standard ai-MBPT scheme^. 

If DFT is used as a reference non interacting system 
Eq. (fl3l) does not hold anymore. Its extension can be 
shown to be 

Cfik ~ £ n k A Ks( nk \ ( r i r i e nk ) 

-Vxc [p\ (?)] I rik) K s, (19) 

where |nk)/<s is the n-th single-particle eigenstate of 
Hks with energy . Note that in Eq. (fTUl) only the 
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£ xc and V xc terms appears as the Hartree terms in £ and 
Vhxc cancels out. 

Eq. 03 reveals the simplicity of the a*-MBPT 
scheme. The accuracy and universality of DFT avoids 
the use of ad-hoc parameters and the prize to pay (at 
least in the electronic case) is to simply subtract from 
the self-energy the xc potential. This simplicity repre¬ 
sents one of the key reasons for the wide-spread use of 
ai-MBPT. 

At this point, one would be tempted to follow the 
same strategy in the electron-phonon case by adding to 
Eq. (1181) . the nuclear Hamiltonian, H n , and the contri¬ 
butions from W e - n that are linear and quadratic in the 
atomic displacements. This is, however, formally not cor¬ 
rect. Indeed, when the nuclei are allowed to be displaced 
from their equilibrium configuration, the DFT (or more 
directly, DFPT) will be expanded in a Taylor series, 

Hks (R) ~ Hks (R) + AH K s (R), (20) 

with, however, 

AH ks (R) ± H (R) - H (R) . (21) 

This is due to the fact that, in AHks (R), the electron- 
phonon interaction terms are statically screened by the 
electronic dielectric function while in the Taylor expan¬ 
sion Eq. (0 they are bare, unscreened. In other words, 
the ground state density p present in Eq. OH actually de¬ 
pends implicitly on the nuclear coordinates. In practice, 
this means that the effect of Vhxc [p] does not appear 
only as an additive term in the Dyson equation but it 
screens the interaction potentials W e - n and W n — n . 

An additional problem, partially connected to the po¬ 
tential double counting of correlation when DFT is used 
as the reference Hamiltonian, is due to the fact that 
most of the electron-phonon theory has been devised in 
the jellium model where the nuclei appear only as static 
and frozen positive charges. As a consequence, strong 
approximations on the perturbative expansion are used 
in textbooks. This is inconsistent with the microscopic 
description of the nuclear lattice and indeed represents 
the most critical shortcoming of the commonly applied 
approaches. In all the aforementioned applications of 
the ai-MBPT schemes (AHC and beyond) W n - n is ne¬ 
glected and W e - n is screened by hand directly in the 
initial Hamiltonian. 

From these simple arguments we can argue that, al¬ 
though az-MBPT is a well-established scheme^ its ex¬ 
tension to the electron-phonon problem is still far from 
being formally defined. 

We would like to apply the MBPT technique to 
the perturbative expansion of Eq. 0 where the non¬ 
interacting Kohn-Sham Hamiltonian and its derivatives 
as calculated by DFT and DFPT are used. 

We propose to apply the standard diagrammatic 
MBPT on the total bare Hamiltonian, given by Eqs. 0, 
explicitly taking into account all interaction terms. We 


will then examine the properties of the A Hks operator 
in order to draw a clear and formal comparison between 
ai-MBPT and DFPT. 

In this way, we aim at creating a consistent framework 
where the role played by screening and self-consistency 
is clearly evidenced even when all the components of 
the electron-nucleus and nucleus-nucleus interactions are 
taken into account. 


III. THE REFERENCE, 

INDEPENDENT-PARTICLE HAMILTONIAN 

As it emerges from the discussion of the previous sec¬ 
tion, the choice of the non-interacting Hamiltonian to be 
used as a reference for the perturbative expansion rep¬ 
resents the connection with DFT and thus provides the 
ab-initio basis for the entire theoretical derivation. This 
is particularly important for the definition of the phonon 
modes. Therefore, we start by introducing a splitting of 
the total Hamiltonian in an independent particle part (for 
independent electrons as well as independent phonons) 
plus interaction terms, 

H (R) = H 0 (R) + AH (R). (22) 

Eq. (1^1) is more suited than Eq. 0 to the MB treatment. 
The reference independent-particle Hamiltonian is 

H 0 (R) = T e + T n + W e _ n (R) 

+ W n . n (R) + AW^J n (R), (23) 

where W e - n and W n - n are evaluated at the equilib¬ 
rium geometry. We have introduced a reference nucleus- 
nucleus interaction Alb"{, a second-order contribution 
within the Taylor expansion in the nuclear displacements, 
that provides the reference phonon modes to be used in 
the diagrammatic expansion. This AW’™{ can be de¬ 
fined from DFPT 

AWZ f n (R) = 

\ E d l^Rc EBO {H)AR lsa AR Vs ,p, (24) 

l s a, V s' fi 

with E bo the Born-Oppenheimer total energy of the sys¬ 
tem calculated within DFT. By construction, the phonon 
frequencies and eigenvectors will be equal to those of the 
Born-Oppenheimer approximation based on the corre¬ 
lated electronic Hamiltonian. Therefore, Eq. (E3l) defines 
an independent-particle Hamiltonian beyond the equi¬ 
librium geometry. The remaining interaction part, up to 
second order in nuclear displacements, is 

AH (R) = W e _e + AW e -n (R) 

+ A W n - n (R) - AW r n e J n (R). (25) 
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where 

A W e _ n (R) = We_„ (R) - W e . n (R) , (26) 

and 

A W n - n (R) = W n - n (R) - W n - n (R) • (27) 

At this point we can introduce the eigenstates of the 
nuclear harmonic oscillations of Hq, written in terms of 
the canonical transformation 

A Ri sa = ^2 (2 NM s ujq X y 1/2 r] a (qA|s) 

qA 

xe i q R is (St_ qA+ 5 qA ) ) (2g) 

where (q, A) is a generic DFPT phonon mode with mo¬ 
mentum q, energy branch A and energy w q A- N is the 
number of q points in the whole Brillouin Zone (BZ). We 
assume the q-point grid to be uniform so that we have 
also N k points for the single particle representation. M s 
is the nuclear mass, rj a (qA|s) is the polarization vector 
of the atom s in the unit cell l in the Cartesian direction 
a, while 6 q a and &q A , respectively, are the annihilation 
and creation operators of the phonon mode (q, A), re¬ 
spectively. 

We now introduce a second quantization formulation 
for the electrons. If </>„k (r) is the H 0 electronic eigen- 


function, we introduce the field operator 


?( r ) = ^E^( r ), 

(29) 

with 


*0k (r) — ^ ^ 0nk (l*) 

(30) 


n 


with c nk the annihilation operator of an electron. By us¬ 
ing field operators, the independent-particle Hamiltonian 
can be written in a second quantization form: 

Ho (R) = El EnkC^kCnk + El ^qA f^qA^qA + ^ J ) (31) 

nk qA ' ' 

where the first term corresponds to T e + (R) and 

the second one to T n (R) + AW"(, (R) . The Born- 
Oppenheimer energy of the ground state at equilibrium 
geometry has been redefined to be the zero of energy 
(hence e.g. W n - n ( R) disappears from this expression). 

The introduction of a reference nucleus-nucleus poten¬ 
tial A WZ2 (R) in Hq (R) was the crucial step to be able 
to map it with its second quantization form. Indeed it is 
well known that phonon dynamics is actually decoupled 
from the electronic one, as shown in many reference s 48 ' 52 . 
Nonetheless, the phonon dynamics can describe accu¬ 
rately the vibrational properties of a real system only if it 
feels the electronic screening. Such screening is accounted 
for in DFPT by the reference potential A W^l n (R) op¬ 
erator. 


IV. THE ELECTRON-PHONON INTERACTION 
TERMS 

Thanks to the definition of the reference part of the 
Hamiltonian, Eq. (1551) , we have that the final splitting in 
independent-particle and interaction terms easily follows 
from Eq. 0. The first two orders of the Taylor expansion 
of H are: 

AH (R) = w(-e + A H {1) (R) + A H {2) (R), (32) 

with 

AR (1) (R) = 

V, d Rlsa [W e _„ (R) + W„_ n (R)]A R lsa , (33) 

z — J lsot 

and 

A H {2) (R) = 

\ E W + Wn-n (R)] 

lsot,l's'fi 

X AR lsa AR Vs , 0 - A W r n e J n (R) ■ (34) 

By using Eq. (E5|) we can manipulate Eq. (1351) and (1551) 
in order to introduce the electron-phonon interaction in 
the basis of the phonon modes. We analyze separately 
the first and second order terms. 

The first order can be written by using Eqs. HMD ® 
as: 



with 


£ q A (r) = d (q A)W e _„ (r, R). (36) 

The function 9( q A )W e - n (r, R) represents the derivative 
of the electron-phonon interaction along the phonon 
mode (q, A). The definitions of the d(q\) operator and 
the £ q A (r) function are given in the appendix [0 Note 
that, in Eq. (l35l) . the real-space integral is performed 
in the unit cell f 0 and not in the whole crystal. This 
is because the sum running on all unit cell replicas has 
been used to impose the momentum conservation at each 
vertex of the interaction terms (see Eq. (1A13I) 1. 

A similar derivation can be done for S q A, which rep¬ 
resents the first-order derivative of the nucleus-nucleus 
potential 

SqA = 0 (qA )W n _ n (R). (37) 

The second-order terms can be worked out in a similar 
way leading to the final form for their contribution to the 
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electron-phonon interaction Hamiltonian: 

A H {2) (R) = 

V / dripl (r) 0 qA , q , v (r) V’k-q-q' (r) 

. k Jo \ 

+ ©qA,q'A' I (ftqX + ^qA^ (^-q'A' + ^q'A'^j • (38) 



In Eq. (l38l) we have introduced the functions 9 and 0 
whose definition is quite similar to Eq. (|3fi|) and Eq. (IdTl) 

0qA,q'A' (r) = p^(qA)(q'A')^e-n ( r i R) > (39) 


QqA.q'A' — n^(qA)(q'A ')W n —n (R) 


-ACi(R) 


(qA)(q'A') 


(40) 


with a w r n e l n (r; 

potential written 
Eq. (l28l) into Eq. 


the reference nucleus-nucleus 

(qA)(q'A') 

potential written in the phonon modes basis by plugging 


The explicit expressions for 9 and for 
^fqA)(q'A') Wn-n (R) are given in Appendix lAl Eqs. (1A18I) 
and (IA19I) . 

Finally, the total Hamiltonian H (R), up to second 
order in the nuclear displacements, can be written as 


r(r) = £ ^nkC nk Cnk + £ f ^qA^qA + 0 + ^ J dridr 2 ft (iq) ft (r 2 ) v (iq - r 2 ) i> (r 2 ) $ (n) + 

nk qA ' ' 


+£ 

qA 


dr^l (r) £ qA (r) ^_ q (r) + £ qA 


(ft qA + K>) + £ [ dr ftk ( r ) QqA.q'A' (r) ^k-q-q' W 

qA.q'A' l/° 


+ 0 


qA,q'A' 


(^-qA + ^a) (' 


&-q'A' + VA' 


)■ (41) 


The diagrammatic transposition of the four electron 
phonon interaction terms in Eq. (I4T1) is given in Fig. ©• 
We see that we have two terms with electronic legs (• 
and ■). Those give direct contributions to the electronic 
propagator. In addition, there are two purely nuclear 
contributions ({££) and □) that do not contribute directly 
to the electron propagator but can be combined with 
the two electronic interactions and still contribute to the 
electronic self-energy, as it will be clear in the follow¬ 
ing. Those are commonly neglected in textbook theories 
of the electron-phonon interaction. But from Eq. (l37l) 
and Eq. 03 we see that there is no reason, a priori, to 
assume that both S qA and 0 qA are zero. 


Gf (t) 

CD a) 




D ( °x it) 

Gib) 




qA 


® = -qA 

Qd) 


The S q A and 0 qA interaction terms do not have elec¬ 
tronic legs because they arise from the purely nuclear 
potential ( W n - n (R)). Nevertheless, at it is evident 
from the above discussion, they can exchange momen¬ 
tum with the electronic subsytem. Energy, instead, is 
not exchanged as the nuclear potential is a static func¬ 
tion. The momentum exchange reflects the change in 
the nuclear-nuclear potential induced by a nuclear dis¬ 
placement. This term is neglected in the jellium model 
because, as it will be clear in the following, the only al¬ 
lowed modes are acoustic excitations for which the zero 
frequency limit corresponds to the zero momentum limit. 
This contribution vanishes as explained in Sec lV A fl 



FIG. 1. Diagrammatic representations of the electron and phonon 
propagators (diagrams (a) and (b)) and of the first (£, diagram (c) 
and S, diagram (d)) and second (0, diagram (e) and ©, diagram 
(f)) order interaction terms in the Taylor expansion of H in powers 
of the nuclear displacements. All interaction terms are written in 
the basis of the phonon displacements. More definitions can be 
found in the text. 























Thus, any coherent and accurate framework that aims 
at providing a consistent way of introducing screening 
and correlation in the perturbative expansion of Eq. m 
will also have to answer to the key question about the 
role played by the nucleus-nucleus interaction. 


V. THE PERTURBATIVE EXPANSION 


non-interacting phonon propagator defined similarly to 
Eq. (|42]> . with bosonic phonon operators 

replacing the electronic ones c: 


D™ (t) = - [1 + n Ka)] (t) + e“« xt 6 (-t) 


- n (w qA ) e^ xt 9 (t) + e u ’* xt 9 (-t) , (47) 


Now that the total Hamiltonian has been split in the 
bare Hamiltonian, Eq. (ED), and in the interaction con¬ 
tributions, W e _ e , AH^ 1 ) (R) (Eq. (l35l) l. and A(R) 
(Eq. (1551) 1. it is possible to perform a standard diagram¬ 
matic analysis. In the following, we will work in the finite 
temperature regime where the electronic Green’s function 
is defined^ as 


G (1, 2) = -Tr{p (0) T t (1) ft (2)] }. (42) 


In Eq. (1551) f3 = ( ksT ) _1 and T is the temperature. T t 
represents the time ordering product, the c„k (t) opera¬ 
tor is written in the Heisenberg representation, p (0) = 
e -p(H-nN) /j , r | e -^(-^-A‘lv)|, w ^h p the chemical po¬ 
tential and N the total electronic number operator. We 
have introduced global variables to represent space and 
time components 1 = (ri,ti). The average spanned by 
the trace operator runs over all possible interacting states 
weighted by the density operator p. 

By using Eq. (1291) , we can expand the Green’s function 
in the electronic basis defined by the reference Hamilto¬ 
nian: 


G(l,2) = 

Y ^«k (ri) 4 >* n , k (r 2 ) G n „' k (ti - i 2 ), (43) 

nn' k 


with 


G„ 


(t) = - Tr[p (/3) T t [c nk (t) 4 k (0)] }• (44) 


The electronic Green’s function can also be expressed in 
a matrix representation as 


[Gk (1, 2 )\ nn , — (/> n k ( r i) <t>n> k ( r 2 ) G nn ' k (ti - t 2 ) ■ (45) 

We will later use the same generic representation of 
Eqs. (l43l) . (l44l) and (l45l) for the self-energy operator. 

The non-interacting electronic Green’s function 
G^k (t) is diagonal in the band index and it reduces 
to a simple exponential^ 


G^(t) = S nn ,e-^-^ 

x [/ (e nk ) 9 (-t) - (1 - / (e„ k )) 9 (f)], (46) 

with / (e) = (eP( e ~G + l) 1 the Fermi-Dirac distri¬ 
bution function. A similar expression holds for the 


with n (e) = (e^ e — l) 1 the Bose-Einstein distribution 
function. Thanks to the standard many-body approach 
and perturbative expansion, it is possible to rewrite 
Gk (t) in terms of G^ 0 ' (t) and the electronic self-energy 
operator Ek by means of the Dyson equation^ 

G k (1,2) = G™ (1,2 ) + j d3d4G^ 0) (1,3) 

x E k (3,4)G k (4,2), (48) 

where the Ek matrix has been introduced following a 
definition similar to Eq. (|45|) . Eq. (|48|) is written in dia¬ 
grammatic form in Fig.EJa). 

Note that in Eq. (|48|) the time variables 1 3 and £4 runs 
in the range [ 0 , f3\. 

In the following subsections we will write E using ap¬ 
proximations with an increasing level of correlation, self 
consistency and screening in order to investigate the ef¬ 
fect of the different electron-phonon interaction terms. 
The solution of the Dyson equation corresponds to an in¬ 
finite series in terms of the Green’s function and the self¬ 
energy. Consistently with the harmonic approximation 
(the expansion with respect to the nuclear displacements 
is limited to the second power in Eqs. 0 and ED) , we 
will work up to second order with respect to A(R), 
and only to first order in A(R). Higher orders of 
nuclear displacements might appear as a consequence of 
self-consistency or screening (in the Dyson equation), but 
we will consider them to be negligible or to have no im¬ 
pact, consistently with our choice of the harmonic ap¬ 
proximation. By contrast, for the electron-electron in¬ 
teraction, there will be no such approximation: higher- 
order powers of the electron-electron interaction will be 
significant. 


A. The Hartree, Debye—Waller and tad—poles 
self—energies 

We analyze first the electronic self-energy that is ob¬ 
tained by considering only one interaction node attached 
to the electronic Green’s function, and select the low¬ 
est non-vanishing diagrams. This self-energy can be ob¬ 
tained mathematically from the Feynmann diagrams in 
figure [2j b, following the diagrammatic rules of Ref. Efll. for 
example. It is composed of four terms: the Hartree T , H , 
Debye-Waller Y , DW , nucleus-nucleus E" _ra and electron- 
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nucleus E e n self-energy: 


£*(l) = y d2v(l,2)G(2,2+), (49) 

Z DW (1) = - E (1) (O'). (50) 

q,A 


E n " n (1) = lim 

q—^0 




(51) 


where u(l,2) = 0 (H — ^ 2 ) v (ri — r 2 ). The electron- 
phonon induced tad-pole contribution is 


-TP, 


r d2^(2)G (2,2+) (52) 


The other diagram (Fan diagram) appearing at second 
order in AH^ (R) will be analyzed in the next sub¬ 
section, Sec. IV B1 together with the Fock diagram from 
W e - e . They both have two interactions nodes attached 
to the electronic Green’s function. 

The self-energy E p is the usual Hartree contribution 
which is a tad-pole diagram made of a Coulomb inter¬ 
action that connects the incoming electronic propagator 
with another, closed, electronic loop. Also the E TPe_n is 
a tad-pole diagram, but, in contrast to the Hartree term, 
the interaction is not electronic but phonon mediated. 



FIG. 2. Dyson equation written in terms of diagrams (frame (a)). 
In frame (b), instead, the lowest-order electronic self-energy in the 
electron-electron and electron-phonon interaction is shown. We 
see the usual electronic Hartree contribution (first diagram) 

the well-known Debye—Waller term (second diagram) 11 , and 

TP _ 

the electron-phonon induced tad-pole (fourth diagram) e n . 
In addition we see the appearance of a new diagram due to the 
derivative of the nucleus-nucleus potential (third diagram), n . 
This diagram plays a crucial role in balancing the contribution from 

TP _ 

£ k e n that, indeed, is not zero in general. 

The nucleus-nucleus self-energy E n_n is a new term 
that has never been discussed before in the literature. 
It comes from the merging of a nucleus-nucleus and 


electron-nucleus interaction. It acquires an electronic 
character thanks to the contraction with the electronic 
propagator embodied in the electron-nucleus interac¬ 
tion (•). 

Finally 'E DW is the well-known Debye-Waller (DW) 
self-energy. It represents the lowest (first) order elec¬ 
tronic self-energy in the second-order derivative of the 
Hamiltonian. 

The total self-energy is local in time and space and 
therefore the Dyson equation of Eq. (H51) reduces to 

G k (1,2) = G< C 0) (1, 2) 

+ J d3G^ (1,3) £ k (3) G k (3,2), (53) 

with 

£ k (l) = 

(1) + X° w (1) + sr n (1) + S k Pe ' n (1) • (54) 

The different contributions shown in Eq. (l54l) have the 
following properties: 

(a) All self-energy contributions are bare. No screen¬ 
ing is present. This is not what should be ap¬ 
plied for practical calculations because, as it will 
be clear in the following, bare potentials lead also 
to unphysical properties. Moreover from DFPT, 
we know that T, DW is screened. In the original 
work of AHC this screening was introduced in a 
semi-empirical manner while in the more advanced 
approach based on DFP T 16 : 57 the electron-nucleus 
interaction is screened in the self-consistent solu¬ 
tion of the Kohn-Sham equation. It is clear, how¬ 
ever, that from a rigorous MB approach this screen¬ 
ing is not present in the original Hamiltonian and 
must be build-up by the electronic correlations. 
How does this screening emerge from a MB per¬ 
spective? 

(b) The nucleus-nucleus self-energy is a new con¬ 

tribution that is not present in any treatment of the 
electron-phonon interaction where the nuclear den¬ 
sity is approximated with an homogeneous charge 
density. In this work, the nuclear coordinates are 
instead coherently taken into account. This is an 
essential step to bridge the MBPT and DFT ap¬ 
proaches. 

(c) In the standard approach to the electron-phonon 
interaction the electron-nucleus self-energy E TPe “ n 
is commonly neglected. However the arguments 
that motivate this approximation^ are based on 
two specific approximations: (i) the nuclear inter¬ 
action is dressed and (ii) there are only acoustic 
modes. In general, however, any system has both 
acoustic and optical modes and the S TPe “ n self¬ 
energy is not vanishing. What is its role and is it 
justified to neglect it? 
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In order to answer those questions we proceed with a de¬ 
tailed analysis of the two series of diagrams connected 
with the dressing of the tad-pole and of the Debye- 
Waller terms. 


1. The electron-phonon induced tad-pole diagram and the 
nucleus-nucleus interaction contribution 

The sum of the E n_n and £ TPe -" is 


This sum would be zero if the £ q A (1) prefactor or if 
the expression between brackets vanishes. However the 
derivative of the bare ionic potential, when q —> 0, di¬ 
verges like |q| _1 . Thus Ea. (l55l) is, actually, divergent. In 
the Jellium model the screening^ of the electron-nucleus 
potential regularizes this divergence and the dressed e-n 
interaction vanishes when q —> 0. This is the standard 
motivation used to neglect the contribution coming from 
the integral appering in the r.h.s of Eg. (1551) . The term 
due to S, instead, has been never considered before. 


S n-n ^ + S TP e _„ 


(!) = lim 

q— 



^a(1)D®(0-) 


X 



d2Qx (2)G(2,2+) 


(55) 


J 


We focus our attention, instead, on the sum of the 
terms between brackets appearing in Eq. (|55l) for small 
but not vanishing values of q. If we work it out we can 
rewrite it in a more clear way. We notice that G (2, 2 + ) = 
p(2) and, from Eq. (1A14I) and Eq. (1A18I) . we have that 


^qA 


dr 2 q A (r 2 , R) G (2, 2 + ) 


r] a (qA|s) e ,q R ' 

h V2 NM s ^ x 


dr 2 W e - n (r 2 , R) p (r 2 ) + W n -„ (R) 


(56) 


Before proceeding in the evaluation of Eq. (156)) we notice 
that from the Dyson equation (Eq. C51) l it follows that 

G k » G k 0) + AG^“ e) + O (A R 2 ) , (57) 

with AG* e_e ) the change in the Green’s function due 
to el-el correlation effects. Eq. m implies that p « 
p(e-e) _|_ Q w ith p( e_e ) the exact charge of the 

system with atoms frozen at the equilibrium positions. 
Therefore, as in Eq. (l56l) d-R. a W e - n (r,R) oc A R, it fol¬ 
lows that, within the harmonic approximation, we can 
safely assume p ~ p ( e_e ). 

It follows than, that within the harmonic approxima¬ 
tion, the overlined quantity in Eq. (1561) is minus the force 
F s acting on the nucleus located at R 


This is true for the exact self-energy but it is not true 
for any approximation of the self-energy unless the Born- 
Oppenheimer energy of the system is calculated accord¬ 
ingly by means of MBPT (for example by using the 
Luttinger-Ward expressions^). 


The condition represented by Eq. (15(1 reveals the cru¬ 
cial role played by the nucleus-nucleus interaction. It is 
only thanks to the coherent inclusion of electron-nucleus 
and nucleus-nucleus contributions that the theoretical 
framework can lead to the justification of the AHC ap¬ 
proach or to the textbooks results. A formal condition for 
the tad-pole diagram to be zero can therefore be defined. 


F/ s = — <9 r. 


dr 2 W e -n (r 2 , R) P {e e) (r 2 ) 


+ W„_„( R) _. (58) 

-I R—R 


By using the Hellmann-Feynman theorem it follows that 

Ftfs = -drUTol H bo (R) |^o) r=S , (59) 


with |4/o) the exact electronic ground state of the total 
frozen Hamiltonian, H (R). 

Thus we can draw the following conclusion: if the self¬ 
energy is chosen in such a way that the nuclear positions 
and density correspond to the exact electronic ground 
state then it follows 

E TPe-„ (-Q + S n-n ^ = 0 _ ( 60 ) 


As an additional approximation we notice that one of 
the most widely approximation used in the litterature is 
to treat correlation effects non self-consistently. In prat- 
ice this means to use the Dyson equation to renormalize 
the single particle energies but not the wave-functions. 
As a consequence, within this approximation, the charge 
density is assumed to be well described by the one calcu¬ 
lated within DFT. This approximation has an important 
and usefull consequence. Eq. m would impose to use 
as ionic positions (R) the ones calculated with a level of 
correlation coherent with the one introduced in p^ e ~ e \ 
As this is a hardly (if not impossible) to do in practice 
the use of the DFT charge allows to approximate both 

p ~ o) and R r; R ( °^ in Eq. (1551) . In this case R® are 
the DFT equilibrium atomic positions that are a simple 
by-product of any DFT calculation. 
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2. Self-consistent diagrams: dressing of the internal 
Green’s functions and of the bare interactions 


Before starting the analysis of the diagrammatic struc¬ 
ture of the self-energy we notice that at any order of the 
diagrammatic expansion we can clearly distinguish be¬ 
tween diagrams that dress the internal electronic propa¬ 
gators and the interaction. 



0 0 
0 + 0 - 0 - 
O “ w 


A clear definition of these two families of diagrams can 
be done by using the simple Hartree approximation for 
the self-energy. Its diagrammatic expansion is shown in 
Fig.© and we notice that at third order, two different 
diagrams appear ((b) and (c)). In the case of diagram 
(b) the interaction builds up a series of bubbles that de¬ 
scribes electron-hole pair excitations. These bubbles, 
when summed to all orders, reduce to the well-known 
Random Phase Approximation for the response function 
as it will be shown explicitly in Sec. (IV A 31) . 


FIG. 3. Diagrammatic expansion of the self-consistent Hartree 
self—energy. Diagram (a) is the non self-consistent contribution 
corresponding to the bare electronic charge. Diagram (b) is com¬ 
posed of bare bubble diagrams. These diagrams belongs to the 
family of diagrams that dress the interaction. Diagram (c), in¬ 
stead, represents a dressing of the internal electronic propagator. 
All diagrams of this kind can be, within the quasi-particle approx¬ 
imations, summed in a definition of a new independent particle 
Hamiltonian, as explained in the text. 


The diagram (c), instead, represents a bare self-energy 
insertion in an internal Green’s function. Any other dia¬ 
gram where an internal propagator is dressed belongs to 
this family. The effect of these diagrams is to renormalize 
the single particle states. This can be easily visualized in 
the quasi-particle approximation where all internal prop¬ 
agator self interaction contributions can be summed in 
the definition of a new set of independent particle ener¬ 
gies, {e„ k }- 

In this work, we are interested in building up a scheme 
to concretely compare the Many-Body and DFT schemes 
as far as the e-n interaction is concerned. In order 
to greatly simplify the analysis, we will focus on the 
first series of diagrams disregarding all diagrams that 
correspond to a dressing of internal electronic propaga¬ 
tors. In the self-consistent Hartree case this amounts 
to neglect the diagram (c) and, in practice, this means 
that the screening of the interaction is described by os¬ 
cillations (described by bubble diagrams) of the bare 
charge. This approximation is commonly used in the 
a?-MBPT scheme and can be written analytically as 


G vaternal _ qM ( 61 ) 


3. Screening of the second-order electron-phonon 
interaction and of the Debye-Waller diagram 

As mentioned earlier, one important aspect that must 
be included in the perturbative analysis in order to bridge 
it with the DFPT formalism is the screening of the 
electron-phonon interaction terms. How does screening 
build up ? This question could appear easy to answer as 
the series of diagrams that screen the lowest order (£) 
interaction is, indeed, well known in the litterature. But, 
what about the second order interaction, S ? 

In the original AHC work, the DW self-energy is writ¬ 
ten, from the beginning, in terms of a statically screened 
W e - n interaction. However this screening cannot be in¬ 
troduced directly in the Hamiltonian. It must appear as 
a result of the diagrammatic expansion. 

In addition, from the discussion of the previous section, 
it follows that E TPe - n (l) ~\~ E” _n (i) j s no t zer0 f or an y 
self-energy that does not reproduce the exact reference 
density and the exact corresponding nuclear positions. 
As a matter of fact this is a condition hard to fulfill in any 
practical implementation as it is computationally very 
difficult to find the nuclear positions corresponding to a 
specific level of approximation for E. 

Even if the condition given by Eq. m can be used as 
a simple approximation it is instructive, for the moment, 
to keep the two self-energies in our derivation in order 
to see their effect on the definition of the screening. 

WvVvV " . + <0000000000 


. .nr FIG. 4. The total bare interaction W„ (1, 2). 

When applied m the diagramatic context we will refer to 

Eq. (1611) as the linearization procedure. The reason of the 

name will be clear in the next section. At the lowest order of the perturbative expansion the 















12 


self-energy shown in Diagram [2jb can be explicitly writ¬ 
ten as 

Sk (1) = E hm | J d2W qA (1, 2) G k (2, 2+) 

+ ^a(i)^S(o-)s; a 

-EW(1 )<!( 0 -)}, (62) 

q ) 

where we have introduced a total bare electron-electron 


interaction W qA (1, 2) (see Fig. (j4])) defined as 

W qA (l,2) = t> q (l,2)J(ii-i 2 ) 

+ ^CqA(l)^° A ) (tl-f2)eqA(2), (63) 

A 

with c q (l,2) the periodic q-component of the Fourier 
transformation of the bare Coulomb interaction 

v (1) 2) = / -^3e^- r ^ q (l,2), (64) 

Jbz (27t) 

with the integral restricted to the Brillouin Zone (BZ) 
only. 

The total bare electron-electron interaction W qA gets 
dressed (i.e. screened) by self-consistency when the 
Green’s function (Eq. (15H10 is inserted into the self 
energy (Eq. (l62l) ). We can also sum over k as the Flartree 
self-energy is a local function and get 


s < 1 > = vl>< 1) = £{ f d2 lim [W q A ( 1 , 2 )] x [g ( 0 ) ( 2 , 2 +)+ f d3G (0) (2,3) £ (3) G (3, 2)1 

k a ^ 

CqA ( 1 ) ( 0 -) S* A j - £ dqA.- q A ( 1 ) D <°> ( 0 -) ( 65 ) 


+ lim 

q —>0 . 


The WqA (1, 2) G k (2, 2 + ) term appearing in Ea. (l62l) is Now we group all S operators to the left hand side of 
written in terms of Feynman diagrams in Fig.([5]a). the equation. We have that 

I 


d3E (3) \s (1,3) - V lim [WqA (1,2)] G<°> (2,3) G (3,2)' 

z ' q —>0 


E { I d2 [W ^ A (1 > G(0) ( 2 ’ 2+ ) + qTo [^ A W D S ( 0_ ) s qA 


E^qA.-qA (66) 

n ' 


and the corresponding diagrams are shown on Fig.([5]b). 


Now, Eq. (16611 is not linear in the sense that the right- 
hand side depends on the dressed G because of the per¬ 
turbative expansion of the inverse of the square bracket 
quantity appearing in the left-hand side. By using the 
discussion of Sec IV A 2l we observe that all dressed G’s 
are internal Green’s functions. This can be deduced also 
by the expansion of the diagrammatic fraction appearing 
in Fig.([5jc). Then we can apply the linearization pro¬ 
cedure , Eq. (16T1) to approximate G(3,2) with (3,2) 
in the Eq. d^Sl) . This allows to define the single-particle 


response function 



Then, we define the dielectric matrix in the Time- 
Dependent Hartree approximation (see later) as 

C(l, 2 ) = 

!,2)-E/ ^3W q A (1,3) x^ 0) (3,2). (68) 

By using Eq. (l66l) into Eq. (l67l) the screening of the 
bare potential W and electron-phonon terms appears so 
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that 


S (!) = E 


d2 lim 

q—>0 . 


+ lim 

q—^0 


Wq\ (1, 2)J (2,2+) 

[^(1)^(0-)^' 


-E^w(l)^S(0-)} (69) 

n ' 


with W, £ and 6 the dressed counterparts of the bare W, 
£ and 8 functions 



€iA (1) = / d2[e^(l,2)]-^ qA (r 2 ), (70) 


$qA,-qA (1) — 


[ d2 

lim (1,2) 

) 

q-m q 


-'qA,— qA 


(r 2 ) 


(71) 


Wq(l,2) = / d3[e^(l,3)] 1 Wq (3,2). (72) 



(0c) 


FIG. 5. Diagrammatic proof that self-consistency at the Hartree 
level is equivalent to screening, at the time-dependent Hartree 
level. The proof is obtained by using the Dyson equation inside 
the definition of the Hartree self-energy. This can be solved in 
terms of the self—energy itself by a simple Fourier transformation 
because the Hartree self-energy is local in time. The mathematical 
transposition of this proof is discussed in the text. 


Eq. (1M1) and Eqs. (1701) . (1711) . (1771) represent an important 
result of the present work. Indeed, they show that self 
consistency screens the interaction lines of all diagrams, 
including the Debye-Waller one. This result will be cru¬ 
cial in discussing how to include higher-order diagrams 
avoiding double-counting problems. 

Indeed, there are two well known ways of increas¬ 
ing the order of the perturbative expansion. One is 
to add skeleton diagrams and the other is to use self 
consistency. This second path is extremely important 
as it provides the way for a given self-energy to fulfill 
conserving conditions^. Skeleton (as well as reducible) 
diagrams are known to build-up the screening of the 
electron-electron and electron nucleus interactions. In 
this section we have shown that screening arises also from 
self-consistency. 


The first two diagrams resulting from the expansion 
of the diagrammatic fraction appearing in Fig. d5]c) are 
shown in Fig. [G] The repeated closed loops represent 
the Time-Dependent Hartree (TDH) approximation for 
the response function. The corresponding screening of 
the interaction is known as Random-Phase Approxima¬ 
tion (RPA). The RPA is the most elemental way to in¬ 
troduce screening in a system of correlated electrons. 

Self-consistency dresses the electron-phonon interac¬ 
tion in the Hartree, in the tad-pole and in the DW dia¬ 
grams. As it will be clear in the following, the equation of 
motion for the corresponding screening function changes 
with the level of approximation used in the self-energy. 
Moreover, when ’£ TPe ~ n (1) + E" _rl (1) ^ 0, the screen¬ 
ing is due to the total time-dependent Hartree dielectric 
function that includes the lattice polarization contribu¬ 
tion. This follows from the definition of the zig-zag inter¬ 
action, Fig. © and Eq. (Rxil) . which induces (see Fig. ©, 
for example) scatterings processes where an electron hole 
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pair is annihilated and a phonon propagator is created. 





©b) 


FIG. 6. First two diagrams contributing to the screening of 
the Hartree (upper diagram) and DW (lower diagram) term. The 
screening can be written as the action of a time—dependent Hartree 
screening function (see text). The dielectric function contains 
only the Hartree exchange diagrams as we considered only the 
Hartree term in the electronic self—energy. However, it also in¬ 
cludes phonon-mediated scatterings as a consequence of the fact 
that T 7 / b--n. has been included in the diagrammatic expansion. 
As it will be clear in the Sec. IV Bl the addition of more diagrams 
to the electronic self-energy corresponds to modify the equation of 
motion satisfied by the dielectric function. 


Even today, most of the calculations at the GW level 
are performed by using the Go Wo non self-consistent ver¬ 
sion. However, this section shows that the DW diagram 
gets correctly screened only by solving the Dyson equa¬ 
tion self-consistently. 


B. The Fock and Fan self-energies 


The analysis of the previous paragraph has been re¬ 
stricted at the Hartree level to keep the notation as sim¬ 
ple as possible. In this section, we extend the derivation 
to the Fock approximation by showing how the screen¬ 
ing of the second-order electron-phonon interaction 6 is 
modified by the inclusion of electronic exchange scatter¬ 
ings via the Fock diagram. We will assume for simplicity 
that E e ~ n + E ra_n = 0 and show the changes that this 
approximation induces in the definition of the dielectric 
function. 

The procedure to disentangle the self-consistency from 
the electron propagator appearing in the Fock operator 
can be done entirely using Feynman diagrams. We start 
by using the Dyson equation to rewrite the dressed (thick 
line) electronic propagators in Fig. ©) in terms of the 
bare electronic propagators (thin line) and of the self 
energy. The resulting diagrammatic expression for E is 
showed in Fig. © a). 



FIG. 7. The Dyson equation at the Fock level in the electron- 
electron and electron-phonon interaction. In addition to the 
Hartree and DW contributions we include the Fock diagram. Note 
that the zig-zag interaction lines include the bare electron-phonon 
interaction defined in Fig. 0. Thus S includes the DW and the 
Fan diagrams with bare and un-dressed interactions. 


We can then work out diagram ©a) by isolating a non¬ 
local operator that is evidenced by the square brackets 
in the diagram ([Sib). The self-energy can be isolated to 
the left-hand side of the equation. Then, in the same 
spirit as for the Hartree case, we can invert the equation. 
By following the same procedure used to go from dia¬ 
gram (j5jb) to ©c) we introduce a diagrammatic fraction 
represented by the square bracket in the diagram ©c). 
Formally speaking, this fraction must be interpreted as 
follows. Let us consider the generic expression -—— T ~. 

with D\ and D 2 two generic diagrams with M open in¬ 
teraction lines (in the present case M = 4). Then we 
have that 


1 

i - Di - d 2 


si + £(D I+ o 2 )", 

n 


(73) 


where the c sub-script means that, at each order, we 
consider the totally contracted products of D\ and D 2 hi 
such a way that the resulting diagram has, again, M open 
interaction lines. We can illustrate this procedure by 
applying it to the DW term (last term in the numerator of 
diagram©c). In this case, Eq. (1731) applied to the square 
bracket produces an infinite series of diagrams that are 
closed in the upper part by a fermion line contracted in 
the second-order bare interaction (■). The first three 
diagrams of this series are shown in Fig. ©a). 

The final result is that, like in Sec IV A 3l the Debye- 
Waller diagram is screened by a dielectric function. How¬ 
ever, there are two important differences with respect to 
the Hartree case. First of all, after linearization of the 
Green’s functions appearing in the right-hand side of the 
diagram ©c) by using Eq. ED, we can define a different 
dielectric function than Eq. (168[) : 


e£ hf (1,2) = 5 (1,2) 


d3vq (1,3) Xq r (3,2), (74) 


with Xq r (3,2) the time-dependent Hartree-Fock irre¬ 
ducible response function. The equation that defines x irr 
is represented in diagrammatic form in Fig. ©b) and 
Fig. ©c). 

The corresponding definition of the screened second- 
order electron-phonon interaction is 


6 , 


qA,—qA 


(i) = 


d2 


lim e^ hf (1,2) 


1 “I 


q—>0 


flqA-qA (r 2 ) • (75) 
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FIG. 8. Diagrammatic proof of the equivalence between self- 
consistency and screening when a Fock and a Fan diagrams are 
present (second and third in the r.h.s. of Fig.JTJ). In the frame 
(b) a portion of the equation is isolated and enclosed in square 
brackets. By defining formally the diagrammatic fraction (see text) 
the equation is inverted and the diagram in the square bracket goes 
in the denonimator of frame (c). The perturbative expansion of this 
fraction leads to the screening of Hartree term and of the second- 
order bare electron-phonon interaction 6. Some of the diagrams 
that follows from the expansion are shown in Fig. ©■ 

The first two orders of x are represented by the two 
bubbles appearing in Fig. 0a) while in Fig. (0 only the 
independent-particle bubbles appear (this is, indeed, the 
definition of the RPA). In this case the Fock and Fan dia¬ 
grams induce a first order bubble with the interaction W 
connecting the two fermion propagators. This diagram 
represents the contribution of the electron- hole attrac¬ 
tion and, when summed to all orders, it can explain and 
predict the formation of excitonic states^. Such bound 
electron-hole states are commonly observed in the ab¬ 
sorption spectrum of several materials^. In this case 
the electron- hole attraction is both electron and phonon 
mediated^. 

The second difference is the fact that in this derivation 
of the Hartree-Fock screening, we have assumed that the 
two tad-poles E e-n (1) and E” _n (l) cancel each other. 
Formally speaking this cancellation is never exact. If 
the contribution of these two terms is included then the 
dashed interaction in the denominator of the diagram¬ 
matic equation 0c) would contain a electron-phonon 
contribution and the whole definition of the screening 
function would be affected. This corresponds, for exam¬ 
ple, to the appearance of exchange diagrams mediated by 


phonons. 



©c) 


FIG. 9. The series (a) corresponds to the first three diagrams 
contributing to the screening of the DW term. Similar series of dia¬ 
grams contributing to the self—energy can be obtained by replacing 
the DW bare self-energy with the Hartree and Fock self-energies. 
As we see besides the usual time-dependent Hartree contribution 
to the polarization function (second diagram) there is a irreducible 
diagram where the electron and hole interact via the total screened 
interaction defined in Eq. 1631) . The definition of the irreducible 
time—dependent Fock response function is given in the diagram (b) 
in terms of the four point function L. The final equation for x > rr 
follows, then, from the corresponding equation of motion for L (di¬ 
agram (e)). 


C. Skeleton diagrams, GW approximation and 
self-consistency issues 

In Sec. IV Al and Sec. IV Bl we have noticed that, even 
if we assume that tad-poles diagrams cancel each other, 
i.e. Eq. (IbUl) is satisfied, the screening of the first (£) 
and second (0) electron-phonon interaction potentials in¬ 
duced by self-consistency depends on the kind of approx¬ 
imation used for the electronic self-energy. 

The situation for the other family of diagrams that 
must be considered at each order of the perturbative ex¬ 
pansion is different. Indeed, if we consider skeleton (bare) 
diagrams we have that the £ function is renormalized by 
the purely electronic dielectric function, as explained for 
example in Ref. [48| (via the diagrammatic method) and 
in Ref. [52| (via the equation of motion approach): 

£ qA (1) = Jd2 [e q (1, 2)\J £ qA (r 2 ). (76) 

In this case [e q | £ ] 1 is the electronic dielectric function 
whose irreducible response function part follows directly 
by contracting the vertex function associated with the 
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self-energy. In the case of the well known GW approx¬ 
imation, the dielectric function is calculated within the 
RPA. The purely electronic expression can be obtained 
from Eq. (1681) when the phonon-mediated exchange con¬ 
tribution is neglected and corresponds to approximate 
WV( 1,2) by Vq (ri,r 2 ). 

(D - 

FIG. 10. The Dyson equation at the GW level in the electron- 
electron and electron-phonon interaction. The electron-phonon 
diagram is known as Fan self-energy and its vertex (represented 
by the circled dot) represents a dressed electron-phonon interac¬ 
tion (see Eq. J76}). The wiggled line is a dressed electron-electron 
interaction (see Eq. CD)- The most important aspect of this di¬ 
agram is that, as long as only skeleton diagrams are included, the 
second-order electron-phonon interaction, and consequently the 
DW diagram, is not screened. 

The non self-consistent self-energy is showed in 

Fig. (USD- The circled dot symbol ( ®) represents the 
dressed interaction defined in Eq. (El and, diagrammat- 
ically, in Fig. ED- Similarly, the wiggled line is the 
screened electron-electron interaction 

W( 1,2) = J d3 [e (1,3)| e; ] _1 v (2,3). (77) 

All the equations and definitions connected to the in¬ 
clusion of skeleton diagrams are well -known in the lit¬ 
erature. The original aspect outlined by the derivations 
presented in the previous sections is that, while the first 
order electron-phonon interaction £ appearing in the Fan 
diagram is screened by skeleton diagrams, the second- 
order 6 is screened by self-consistency. 

This deep difference in the procedure that defines the 
kind of screening of the electron-phonon interaction is re¬ 
flected in the different equation that is satisfied by e in, 
for example, Eq. (1771) , Eq. (1751) and Eq. (1701) . Depending 
on the choice of the self-energy, we have phonon me¬ 
diated exchange and/or direct scatterings and electron 
hole attraction diagrams. We will see in the following 
that yet another family of dielectric functions is used 
within the DFPT approach. The physical interpretation 
of these different definitions is discussed in Sec. IVIIBI 



FIG. 11. Diagramatic representation of the dressed electron- 
phonon vertex within the GW approximation for the self-energy. 
In this case the diagrams are bare (skeleton) and sum into an RPA 
dielectric screening of the ionic potential. In this case, even in 
the case where the tad-pole diagrams are non zero, the dielectric 
screening is purely electronic. The standard additional approxima¬ 
tion is to consider a static screening. 


VI. THE DENSITY-FUNCTIONAL-THEORY 

APPROACH TO THE ELECTRON-PHONON 
PROBLEM 

In the previous section we have analyzed the kind of 
diagrams induced by the electron-electron interaction in 
the dressing of the electron nucleus interaction terms. 
We have disclosed the key role played by self-consistency 
and the different level of approximation that arises from 
the perturbative expansion. 

In order to link the DFPT and MBPT approaches we 
start with a short review of the purely DFT-based ap¬ 
proach to the electron-phonon coupling. 

DFT is a self-consistent theory, and DFPT is its exten¬ 
sion to take into account, sclf-consistently, the effect of 
static perturbations (like nuclear displacements). In this 
case, DPFT provides an exact description of phonons 
within the limits of a static and adiabatic approach. 
The phonon frequencies in DFPT are always real and no 
phonon dissipation process is included. If we introduce a 
total (electron-electron plus electron-nucleus) potential 

V scf (R, r) = V Hxc [p] (r) - ]T , (78) 

u l r - R ^l 

where the functional dependence of the electronic density 
on the nuclear positions introduce a direct (via W e - n ) 
and indirect (via Vhx c) dependence on (R) in 14c/- In 
DFPT this complicated dependence links the calculation 
of <9r ls V sc f (R,r) to the solution of a self-consistent set 
of equations. 

By applying the same procedure used to derive Eq. 0. 
a formal Taylor expansion of Hks can be obtained. How¬ 
ever, if the dependence of the density on the nuclear posi¬ 
tions is not taken into account, all terms in the Taylor ex¬ 
pansion are bare. In DFT (or DFPT), screening builds- 
up because of the V sc f indirect dependence on (R). We 
introduce 

£q\ PT ( r ) = d(q\)V sc f (R, r), (79) 

and 

CTTa (r) = 0 (q A),(- q A) V scf (R, r). (80) 

The expression for £, PPPT (r) can be written in terms of 
<9r 3 V ac f (R, r) by following the same procedure outlined 
in appendix [A] The screening of £ within a pure DFPT 
scheme follows from the fact that 

<9r i B V sc f (R,r) = <9 r ls W e -n (R,r) 

+ / (81) 

In order to evaluate Eq. m and create a link with the 
MB approach, we notice that DFPT is based on the linear 
response regim o 16 ' 57 where: 

dR u p(r) = J dr'xDFT (r, r 7 ) <9r 1s W e _„ (R,r'), (82) 
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where the DFT polarizability xdft = g V p is solution 
of the following Dyson equation 


Xdft (r,r') =X { ks ( r > r ') 

+ Jdr"^ { r,.)f H ^y)XDFT(r"y), (83) 


with Xk\ th e independent particle KS response function. 

From the definition of the Hartree and xc potential it 
follows that 


fnxc (r,r') 


SVhxc [p\ (r) 

5p{v') 

v (r - r') + f xc (r, r 7 ), 


(84) 

(85) 


which, finally, yields the well-known expression for the 
derivative of the total DFT self-consistent potential 


dn la V sc f (R, r) = j dr'{<5(r-r') 

+ J dr" f H xc ^5 ^ X.DFT ( r "’ r )} 

dn lB W e - n (R, r'). (86) 

If we now introduce the DFT dielectric function 
[ e UFT ( r , r ')] _1 = J( r -r') 

+ J dr" f H XC ^^*5 ) Xdft (r ,r') , (87) 
we have, finally, that 

dR la V scf (R, r) = J dr' [e DFT (r, r')] 1 

xa Rls lF e _„(R,r'). (88) 

Similarly, the second-order derivative of V sc f can be in¬ 
troduced-^^ as 

^R !s R,/ s /^sc/ (Rj r ) = 

dn ls |y dr' [ e DFT (r, r')] 1 9 r, (( i , W e - n (R, r')j. (89) 

Eqs. (1551) and (1551) must be compared with Eqs. G0, CD 
and d72l) in order to highlight the differences between the 
two formulations and potential similitudes. 

Eq. (1551) can be written in the basis of phonon displace¬ 
ments (q, A) as 

fqA FPT « = J dr' [e FFT (r, r')] _1 £ FFPT (r') • (90) 

This last equation can directly be compared with Eq. GOD 
and the following observations can be made 

(a) In DFPT the £ function (and more generally the 
electron-phonon interaction) is statically screened 


and it does not include the contribution from 
the lattice polarization. In the MBPT, instead, 
the electron-phonon interaction is dynamically 
screened (i.e. the dielectric functions defined in 
Eq. (l68l) and Eq. d74l) are time dependent) and it in¬ 
cludes phonon mediated scatterings (see Eq. (1701) ). 
The static screening of DFPT reflects the fact that 
there are no retardation effects in the theory. Those 
effects are peculiar of the MBPT and, in some 
cases, can lead to important deviations from the 
static limi t 11 ! 59 ' 60 when included in the self-energy. 

(b) The DFT dielectric function e DFT defined in 
Eq. (1571) is a test-electron dielectric function 
whereas in the MBPT the dielectric function that 
screens the bare electron-phonon interaction is a 
test-charge function. The difference between those 
two functions is well described in Ref. (gTI and |62j. 
In the test-charge case, the dielectric function rep¬ 
resents a response to an external particle, while in 
the test-electron case, the charge is itself an elec¬ 
tron. We will discuss in more detail this difference 
from a physical perspective in Sec lVII 51 

(c) A peculiar consequence of the DFPT approach is 

the appearance of non-rigid nuclei contributions 
to the second-order derivative of the electron- 
nucleus interaction potential, ,V sc f (R, r). 

This contribution does not appear in the Many- 
Body derivation carried on in the previous sections. 

Such non-rigid nuclei (recently called in the literature 
non rigid ion contribution^) can be studied by applying 
the derivative with respect to a nucleus displacement on 
the right-hand side of Eq. (1551) and by distinguishing a 
rigid-nuclei (RN) and a non-rigid nuclei (NRN) terms 


^R is Ri/ s / Vscf (R, r ) — t , R, s R i , s , Vscf (R, r ) 


RN 


^Ri a R lla ,Vscf (R, r ) 


NRN 


, (91) 


with 


^Ri.Ri'.'^W (R’ r ) 


RN 


J dr'[e DFT (r,r')] 1 ^ IsRi , a ,W e _„ (R, r'), (92) 


and 


^Rf -jR,/,/ Vscf (R,r) 


NRN 


J dr'd Rl5 [e DFT (r,r')] 1 <9 Ri , a , W e -„ (R, r'). (93) 

If we now rewrite both terms in the DFPT phonon (q, A) 
basis and multiply by a 1/2 pre factor we get 


nDFPT,RN , x _ 
°qA,q'A' PI — 


1 


, -1 


dr' [e DFT (r,r')]~ 

x^A.q'A'VWtR.r'), (94) 
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and 


A. Tad—Poles 


qDFPT,NRN , N 
AjA.q'A' VC — 

\ J dr'^r (r') dqx [e DFT (r, r')] _1 . (95) 

From the local dependence on R of the V sc f it follows 
that, in the RN contribution, a dR ia> R 1/s , appears. In 
the rigid-nuclei approximation (also called the rigid-ion 
approximation) the d RFR , R ’, NRN is neglected. 

The physical interpretation of the NRN contribution 
to the second-order derivative of the self-consistent po¬ 
tential is obscure in the DFPT derivation and seems to 
be more a mathematical separation based on computa¬ 
tional load of the calculations rather than a physically 
motivated choice (see Ref.for a detailed explanation). 
We will discuss this term from a Many-Body perspective 
in Sec. IVII Cl 


VII. MBPT STARTING FROM 
DENSITY-FUNCTIONAL AND 
DENSITY-FUNCTIONAL PERTURBATION 
THEORY 


The main difficulty in performing a diagrammatic ex¬ 
pansion on top of DFT is that it is not possible to write 
the initial Hamiltonian in terms of dressed interactions. 
This means that, even within DFT, the Taylor expansion 
of H (R) is still given by Eq. (1321) with the only differ¬ 
ence that the electron-electron interaction is replaced by 
the mean-field xc-potential. At this point DFPT and 
MBPT follow two different routes in order to describe 
the dressing of the interaction and the definition of the 
electronic self-energies. 

To obtain the same screening as DFT+DFPT from 
a purely many-body point of view, we start by writing 
again the total Hamiltonian as a function of the KS one 

H (R) = H ks (R) + H n (R) + A W r n e J n (R) 

Y A H rerna i n i n g (R) . (96) 

We can determine the value of A H rema i n i ng (R) from 
Eqs. (fT71) . (031) and C5l) 

A H rernain in g (R) = A W n . n (R) - AW^ (R) 

+ A W e -n (R) + W e _e - V Hxc [p\ , (97) 

with ~p the equilibrium DFT density. 

Now the full MBPT machinery described in the pre¬ 
vious sections can be applied to Eq. leading to the 
screening of the electron-phonon interactions. However, 
our aim is to create a link with the DFPT definitions 
of Eq. ([88]) and (l89l) . The main differences between the 
MBPT and DFPT approach are listed in Table Q] and 
explained below. 


By definition the DFT density is the exact one and it 
satisfies Eq. (1581) . Thus, bearing in mind the linearization 
procedure outlined at the end of Sec. IV ATI we can affirm 
that within DFPT, £ e_n (1) -f £™ _ " (1) = 0. As a conse¬ 
quence the dielectric function that screens the electron- 
phonon interaction does not include phonon-mediated ex¬ 
change scatterings. 


B. Screening of the second—order electron—nucleus 
interaction 


In order to discuss how screening builds up, we start 
from the lowest order self-energy in the electron-phonon 
scattering (see Fig.J2]b)). Again, we start from the bare 
Hamiltonian but this time we use the KS system as refer¬ 
ence. The Dyson equation that follows from the Hamil¬ 
tonian of Eq. (1TO1) is 


G (1,2) = G ks (1,2) + J d3G KS (1,3) 

x (3)-V Hxc [p\(3)}G(3,2). (98) 

Moreover, since the Hartree part (Vjj) also appears in 
Vhxc , we have 


[ P ](i)-v H [p](i) = 


d2v (1, 2) p(2)-p(2) . (99) 


In the case of a local self-energy, the difference of densi¬ 
ties can be rewritten in terms of the self-energy by using 
Eq. ((53]) 

P (1) - P (1) = J d2G KS (1,2) £ (2) G (2,1). (100) 

Following Eq. (1551) . we obtain within RPA and after lin¬ 
earization of the Green’s function (G ~ G KS ) that 

e KS ’ RPA (1)2) = £ (lj2 ) J d 3v (1,3) \ KS (3,2), (101) 

with 


\ KS (1, 2) = G ks (1, 2) G ks (2,1). (102) 

The final expression of the second-order electron-phonon 
interaction thus becomes 


^qA,-qA (1) 


d2 


lim e 

q—^0 


KS,RPA 

q 


( 1 , 2 ) 


-1 


X 6»qA,—qA (r 2 ) 5 (f 2 ) , 
with the dielectric function defined by 
[ e KS,RPA^ 2 j j- 1 = ,5(1,2) 

+ j d3v(l,3)x KS ’ RPA (3,2), 


(103) 


(104) 
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TABLE I. Schematic representation of the different treatment within MBPT and DFPT of the most important aspects of the 
electron-nucleus interaction. 



MBPT 

DFPT 

Tad-Pole diagrams sum, E e “ n + E n “ n , 
Fig. J2]b). 

Vanishing only when the atomic positions 
are evaluated consistently with the level 
of correlation included in the self-energy. 

Vanishing. 

Screening of the first order interaction, 
£ ( r )- 

Time-dependent and described by a test- 
charge dielectric function consistent with 
the self-energy. Induced by skeleton 
diagrams. 

Screened by the static test-electron DFT 
dielectric function. Introduced via self- 
consistency. 

Screening of the second order interaction, 

e(r). 

Time-dependent and described by a 
test-charge dielectric function consistent 
with the self-energy. Induced by self 
consistency diagrams. 

Screened by the static test-electron DFT 
dielectric function. Introduced via self- 
consistency. 

Non-rigid nuclei contribution to the sec¬ 
ond order interaction, 6 (r). 

Appears in the static limit from the di¬ 
agrams that describe the dressing of the 
electronic density caused by the e-n in¬ 
teraction. It is, again, caused by self 
consistent diagrams. 

Present and static. Induced by the im¬ 
plicit dependence of the screening on the 
atomic positions. 


and x KS ’ RPA solution of Eq. (1551) with fuxc (r, r') = 
w(r, r'). 

We notice that Eq. (11041) defines a test-charge dielec¬ 
tric function while in DFPT, Eq. (1571) . a test-electron 
dielectric function appears. 

This difference of definition is a straightforward con¬ 
sequence of the fact that the definition of the dielectric 
function is linked to the distinction between the classical 
and the quantum parts of the induced potential. The 
classical part satisfies a Poisson equation whose solution 
is the Hartree potential. The quantum part is treated 
in different ways in MBPT and DFPT. In the MBPT, 
the quantum induced field is represented by the change 
of the correlation self-energy due to a test charge. This 
is described by the vertex function^ that can also be 
used to rewrite the exact self-energy in a closed form. 
In DFT, instead, electronic correlations are included in a 
mean-field manner by means of the exchange-correlation 
V xc potential. It follows that, the variation of V xc mimics 
the variation of the self-energy and, thus, represent the 
quantum part of the induced potential. Therefore, the 
difference between a test-electron and a test-charge is 
that the test-electron includes the total variation of the 
total potential, including V xc . This contribution leads to 
the f xc term appearing in the right-hand side of Eq. © 
and marks the difference between the MBPT and DFPT 
screening. 

As additional proof, we can notice that if the electronic 
self-energy is approximated with a DFT exchange- 
correlation potential, then the many-body vertex turns 
the test-particle into a test-electron dielectric function 
that is consistent with the DFPT definition. This means 
that f xc is taking into account, in a mean-field manner, 
the effect of the MBPT vertex function. 

At this point, we can conclude by observing that if 
local or semi local approximations for V xc are used, the 
difference between a test-charge and a test-electron di¬ 
electric function can be safely neglected. Indeed any local 


or semi-local expression for V xc is regular in the short- 
distance limit whereas v (r, r') diverges. This means that, 
when q —> 0, it follows that Jhxc (r, r 0 ~ v (r, r'). 

More elaborate expressions for V xc that also include 
proper short-distance spatial corrections exist. In this 
case, a more accurate analysis of the effect of the Many- 
Body vertex function on the screening of the electron- 
phonon interaction becomes essential to draw a conclu¬ 
sive parallel between MBPT and DFPT. However, this 
goes beyond the scope of the present work. 

Therefore, as far as local or semi-local approximations 
for V xc are used, we can approximate the DW screened 
interaction with the one evaluated at the DFPT level. 


C. The non-rigid nuclei contribution to the 

Debye—Waller self—energy from a diagrammatic 
perspective 

One of the most undeniable difference between the 
DFT and MBPT scheme is the absence of a diagram¬ 
matic explanation for the non-rigid nuclei contribution 
(Eq. ®) to the second-order derivative of the nuclear 
potential. This terms has been shown to be quite impor¬ 
tant in low-dimensional systems^! but, at first sight, it 
does not appear in the many body theory of the electron- 
phonon interaction. Indeed, in Eq. ED and Eq. ED. the 
derivatives with respect to the atomic positions act only 
on the local bare ionic potential. 

In order to find a diagrammatic perspective of 
^qA q'A' we foll° w again the path of performing a 
diagrammatic expansion of the bare, un-dressed Hamil¬ 
tonian and, at the end, draw links with DFT. In the 
present case we need to to consider a new series of di¬ 
agrams describing the change of the electronic density 
induced by the electron-phonon interaction. Three ex¬ 
amples of this series of diagrams are showed in Fig. ED. 
We can see that those diagrams are of the same order of 
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magnitude as the lowest order self-energies and there is 
a priori no reason to neglect them. 



Clb) 


In order to introduce a simple and clear interpreta¬ 
tion of the contribution due to the 5T, H diagrams, we 
take the static and adiabatic limit of the atomic dis¬ 
placements. This approach will also simplify the con¬ 
nection with the corresponding quantity evaluated with 
the DFPT scheme. 

We start by taking the static limit of ( t 3 — i 4 ) ~ 

S (t 3 - U) [2n (w q A) + 1] and of £(3) « £(r 3 )£(f 3 ). 
Moreover, in this limit, we can treat the atomic dis¬ 
placements as classical and static variables in order to 
approximate AH (R) 

A#(R)«^ e q A (r.,:)n qA , (109) 

qA,i 

with Uq\ the phonon displacements defined as 



me) 


FIG. 12. First three diagrams that dress the Hartree self-energy 
due to the change of the electronic density induced by the electron- 
phonon interaction. The sum of all diagrams of this kind (see 
Eq. (11051) for a more formal definition) reduce, in the static and 
adiabatic limit, to the NRN contribution to the Debye—Waller self¬ 
energy as discussed in the text. 


The dressing of the Hartree self-energy can occur in 
two ways: the first is by dressing the internal electronic 
propagators (an example is the Diagram (fili al). The sec¬ 
ond way is to dress two tad-poles with phonon scatterings 
like in Diagrams (fl2l b) and (fl2l c). All the diagrams of 
this kind can be written as 


Uq X = ^2 ( 2NM sUq\) 1 AR lsa rj a (qA|s). (110) 

Isa 


As we are interested in a specific series of diagrams, we 
can disregard higher-order corrections to AR(R). The 
Eq. (11091) allows to formally define the derivative of the 
Green’s function, dq\G(l,2) evaluated at the equilib¬ 
rium nuclear positions. Indeed, from Eq. (11091) . it follows 
that 


[GM ]" 1 


G (0) (1,2) 


AH(R). (Ill) 


and 


< 9 qA [G ( 1 , 2)]- 1 = 

— dq\AH (R) = —£qA ( r i) S (1,2), (112) 

which, finally, using the identity 

6(1,2) = J d3 [G (1,3)] _1 G (3,2), (113) 


5T, H ( it) = J dr 2 v(r 1 - r 2 )Sp(r 2 ), (105) 

with 

6p (r 2 ) = \p (r 2 ) - p (r 2 )] = ^ 6p p (r 2 ) (106) 

P 

= ^ ^n'p fo) 6pnn' p (r 2 ) ■ (107) 


yields the desired definition 

dq X G(l,2) = J d3G^(l,3)f q A(r 3 )G® (3,2). (114) 

Eq. (|114|) is diagrammatically represented in Fig. m- 
The derivative of the Green’s function splits the projector 
in two with the insertion of a dressed electron-phonon 

first-order interaction ( ® ). 


In the case of the diagram (|12l a) we have that 




(§> 


Sp p (2) = J 2 / dM [^p 0) (2) 3) C q A (3) 


q,A 


x Gp_ q (3,4) q A (4) (4,2) (t 3 - t 4 ), (108) 


FIG. 13. Diagrammatic transposition of the derivative of a Green’s 
function with respect to a specific displacement written in the 
phonon basis (see text, in particular Eq. (11141) '). The derivative 
splits the Green’s function with the insertion of an electron-phonon 
first-order interaction. 


where £ is a first-order interaction that we assume to 

be screened by the very same skeleton diagrams that we By using Eq. (11141) we notice that the first-order inter- 
described in Sec lV Cl actions appearing in the p (r 2 ) — p (r 2 ) can be interpreted 
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as derivatives of the Green’s function. By using Eq. fl 1141) 
we can indeed rewrite the quantity in square brackets of 
Eq. (11081) as 

J d3G™ (2,3K qA (r 3 )G^ q (3,4) = 

d qA G k (2,4). (115) 

We take now the static limit of y4°), which is consistent 
with the static and adiabatic approach taken in this sec¬ 
tion. In this way, we obtain the contribution of the first 
diagram ((T2l al to 5T. H 

5T, H ’ a (ri) = ^2 I rfr 2 dr W ( r i,r 2 ) 

^ qA ■’ 

x 9 qA x (0) (r 2 , r 3 ) f* A (r 3 ) [2 n (w qA ) + 1], (116) 

where the factor 1/2 follows from the definition of 
(see Eq. (1671) ) that yields two equivalent contributions to 
6Z H ’ a . 

If we now perform the same procedure to all bub¬ 
ble diagrams that contribute to time-dependent Hartree 
reducible response function we get the diagram (fl2lb'l 
and (He). The final expression for S'S 11 is 

<5£ ff ( ri ) = ^ / dr 2 d^ [e tdh (r 1 ,r 2 )] _1 
qA J 

x CqA (r 2 ) [2 n (w qA ) + 1]. (117) 

We have demonstrated that, in the static and adia¬ 
batic limit, the series of diagrams due to the dressing of 
the electronic density as a consequence of the electron- 
nucleus interaction reduce to a new contribution to the 
series of second-order diagrams. However, the derivative 
with respect to the atomic displacements is acting on the 
dielectric function. 

The Eq. (11171) allows us to rewrite the DFPT NRN 
potential 9q\qf'x’> NRN in a purely MBPT language. A 
direct comparison between Eq. (11171) and Eq. (l95l) reveals 
that the difference between the MBPT and the DFPT 
approaches lies in the dielectric function (e td?l for the first 
and e DFT for the second). 

D. Screening of the first-order electron—nucleus 
interaction 

If we now consider terms beyond the Hartree approx¬ 
imation, the definition of the dielectric function that 
screens the second-order interaction will change. This 
means that if we add to the Hartree potential a Fock 
operator, a Vh xc potential or a GW self-energy, the di¬ 
electric function will solve a time-dependent Fock, DFT 
or GW equations. 

We can anyway add new correlation terms to the self 
energy by keeping frozen the self-consistency at the DFT 


level. This means that the dielectric function will approx- 
imatively solve a DFT equation that, within the limits of 
the discussion of Sec. IVIIBI will correspond to a simple 
RPA approximation written in the KS basis. 

If we now include the Fan contribution to the total 
self-energy and follow the same procedure already dis¬ 
cussed for the skeleton diagrams (see Sec. IV Cl) , we will 
obtain that the £ function is renormalized by the stan¬ 
dard skeleton polarization diagrams. These, again, will 
be summed in a RPA dielectric function written in terms 
of KS bare Green’s functions. 


E. Self-consistency 

As already mention earlier, the self-consistent dia¬ 
grams play a crucial role in the definition of the second- 
order interaction. 

The screening of the rigid-nuclei second-order interac¬ 
tion can be obtained with a simple self-consistent Hartree 
approximation (that is already embodied in the DFT pro¬ 
cedure). Instead, the non-rigid nuclei term requires to 
include, self-consistently, Fan diagrams and cannot be 
obtained with a simple mean-held approximation. 

This complicates the merging of MBPT and DFT as it 
is now clear that we can approximate the 6 and £ func¬ 
tions with the DFPT counterparts under the condition 
that the Dyson equation is not solved self-consistently. 
Any kind of self-consistency would screen again the 9 
interaction leading to a severe over-screening. 

The self-consistency must be taken from the solution 
of the KS equations and not re-introduced at the MB 
level. This, as explained above, will correspond to the 
inclusion of certain class of diagrams, a practice that is 
well motivated within a many-body approach. 

Finally, if the Dyson equation is not solved self 
consistently, any additional electronic self-energy beyond 
the V xc can be introduced as an additive potential (see 
Eq. (flbl) '). This means that it is safe to introduce non 
self-consistent quasi-particle corrections (electron and/or 
phonon mediated). 


F. Higher—order derivatives of the V sc f 

The importance of creating a coherent merging of DFT 
and MBPT stems from the possibility of deriving an ac¬ 
curate and predictive approach to the electron-phonon 
interaction. There is however an entire family of physical 
problems induced by a strong electron-phonon interac¬ 
tion that requires to introduce higher-order self-energy 
diagrams. A crucial difference with respect to the state- 
of-the-art models (like Frohlich or Holstein Hamiltoni¬ 
ans) is that the correct treatment of the nuclear positions 
produce a new category of diagrams of generic order n 
where the n-th order derivative of the self-consistent po- 
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tential appears 


rp{n) 

± Ry—R n 


= d n V scf (R, r) 

{ ’ <9Ri • • • <9R n ’ 


(118) 


where the Debye-Waller diagram is the lowest order ex¬ 
ample. Thanks to the simplicity of DFPT, it is possible 
to define an iterative expression for the matrix elements 


of T, 


(n) 


Ri"-R 7 


(r) 


T, 


(n) 


Rl ■ ■•Rn 


nn ' = ( nk \ T R?-R n ( r ) l n ' k ')- 

kk' 


(119) 


Such matrix can be rewritten in terms of T^ n ^ and of 
the n-th derivative of the KS energy levels 


rp(n) 

J -R 1 —R n 


7 — Jr. 


nn \ s ^mp.nk 
kk' 


7^(1) rp{n— 1) 

± Ri nm R 2 ■ “Rn 

kp 


mn 

P k " 


^mp ^nk 


rp(n—l) 

J R 2 -R n 


7~>(i) 

nm Ri run' 

k P_pk / 


^n'k' ^mp 


+ Pi 


( rp(n—l) 

I 1 R 2 -R, 



( 120 ) 


with 


rp(0) 


mn 

pk' 


— ^n'k'^n'm^k'p. 


( 121 ) 


and 


PiF (Ri... R„) = ~id Rl F (R x ... R„), (122) 


where F is a generic function of the atomic positions. 

More information on Eq. (11201) is provided in ap¬ 
pendix E Such procedure can be used to calculate nu¬ 
merically, in an ab-initio manner, the high-order self 
energies. 


VIII. CONCLUSIONS 

In this work we have studied the electron-phonon 
problem by comparing the standard Many-Body Pertur¬ 
bation Theory (MBPT) with Density Functional The¬ 
ory (DFT) and Density Functional Perturbation Theory 
(DFPT). By analyzing the different diagrams that con¬ 
tribute to the electronic self-energy we have achieved sev¬ 
eral important goals. 

(i) The well-known electron-phonon induced tad-pole 
diagram is not zero in general but can be cancelled by 
a nuclear-nuclear self-energy if the equilibrium nuclear 
positions are coherent with the level of correlation in¬ 
troduced at the single-particle levels. In the case of a 
non self-consistent calculation, the equilibrium positions 
evaluated with a DFT reference Hamiltonian removes the 
sum of the tad-pole diagrams. 

(ii) Self-consistency diagrams dress the second-order 
interaction and the corresponding Debye-Waller self 
energy. This provides the many-body interpretation of 


the screened Debye-Waller self-energy already known in 
the DFPT case. 

(iii) We identify the specific series of diagrams that 
explain the non-rigid nuclei contribution to the Debye- 
Waller self-energy. The existence of this term was known 
only in a purely DFPT (static and adiabatic) approach. 
In the present work, we provide a clear physical inter¬ 
pretation of this term by performing a static limit of the 
MBPT expression to demonstrate that, indeed, it reduces 
to the well-known DFPT result. 

(iv) We have drawn a final series of statements re¬ 
garding the possibility to perform many-body pertur¬ 
bation theory calculations on top of Density Functional 
and Density Functional Perturbation Theory avoiding 
the double counting of diagrams. 

This work represents a firm and formally accurate in¬ 
spection of the two methods while describing the lim¬ 
itations of the static DFPT approach and providing a 
practical way to go beyond by the merging with more 
advanced MBPT methods. 
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Appendix A: The interaction Hamiltonians in the 
phonon displacements representation 


We proceed here to discuss the details of the derivation 
and definition of the functions 

£qA (r) = d(qA) W e -n (r, R) , (Al) 

^qA.q'A' (r) = q^fqA^q'A') W e -n (r, R) j (A2) 

SqA = <9(q \)W n - n (R) , (A3) 


©qA.q'A' — r,^pjA)(q'A ')Wn-n (R) 


- A W?J n (R 


(qA)(q'A') 


(A4) 


where we purposely did not place a nucleus dependence R 
on the left-hand side of those equations. We will indeed 
see that the dependence is lifted by Eq. (IA9D . as the 
derivatives have to be evaluated at R. 

Those functions enter in the different interaction terms 
of the electron-phonon interaction . We start by remem¬ 
bering that 


W e _„( r ,R) = -^Z s <;(r-R, s ), (A5) 

Is 
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and 

W n -n (R) = l V' „ , Z s z a ,v (R u - R Va .). (A6) 

Z Z - * Is,l's' 

It follows that 

d Rlsa W e _ n (r, R) = Z s v (r -R ls ), (A7) 

and 

d Ru a Ri>s>«’ We - n ( r ’ R ) = 

-Z B v™ (r-R ls )5 w 6 SS ', (A8) 

with (r) = d ra v( r) and v^, (r) = dl ar ^v(r). By 
using Eqs. m and (IA8I) the evaluation of the <9( qA ) and 
9^ qA )( q , A ,) is straightforward. It is enough to remember 
that 

5 (qA) E (R) = E (27VM sWqA )- 1/2 

Isa 

x Va (qA|s) e^d Rlaa F (R)| r=r , (A9) 

with F (R) a generic function. The resulting expression 
is therefore independent of R. The c>^ qA )( q ' A ') is then 
obtained by applying twice Eq. (IA9I) . We finally obtain 
that 

£ q A (r) = 

(A10) 

Isa V ZI ^ MsUJ 1 x 

The same machinery can be applied to d qA . q ' A ' (r) ob¬ 
taining: 


the sum in the unit cell corner. It follows that Eq. (IA12I) 
turns into 


a 1 U" r s 


e *(q+P" k ) R i 

L l 


x 


E [ dr < k « ( r - r «) ( r ) e l(q+p “ k) r , 

sa ^ ® 

(A13) 


where u n k is the periodic part of the wavefunction. By 
using the fact that ^ e lPRi = NSp we finally obtain 
that 


(nk|£ qA (r) |rc'p) 


E 


= d P ,k- q (nk|£ qA (r) |n'k 

Zse^-VN , ,, , 

Va (qA|s) 

y2M s w qA 


q) = 


/ * <k (r) vi 1} (r — r s ) « n - k - q 0) 


(A14) 


L^O J 

The same strategy can be used to define the matrix ele¬ 
ments of 6: 


(nk|d qA . q / A / (r) |n'p) = 

<5p,k-q-q'(nk|d qA!q / A / (r) \n'k - q - q') = 

- E fE va ( qA i s ) 7?q ' ( q,A 'i s ) 

w V 4M ^qA 

/ dr < k (r) v™, (r - t s ) u n , k -q- q ' (r) 

Vo 


(A15) 


The definition of S qA is similar and follows directly from 
Eq. (HU) and from the extension of Eqs. (1A71) and Hi 
to the nucleus-nucleus potential: 


y p i(q+q')-R.!s 

<Wv (r) = - E ' 

Isaa' 2NM S y/Wq\U}qi\> 

x r] a (qA|s) r] a , (qA|s) v™, (r - R (s ) . (All) 

Two important quantities follow from the integral 
of Eq. (IA10I) and (IA11I) when multiplied by two 
single-particle wavefunction: (nk|£ qA (r) |n'p) and 

(nk|d qAiq / A / (r) |n'p). In order to evaluate those we have 
to consider, in the first-order derivative case, a term like 


E 


dr C k (r) v R) (r - R is ) <t> n , p (r) e lq ' 


iq-Ris 


, (A12) 


with the r integral performed on the whole crystal. Now 
we observe that i/ 1 ) and v^ depend on the atomic posi¬ 
tions only via their argument. If R/ s = Rj + r s with t s 
the position of the atom s inside the unit cell located at 
R, we can change variable from r to r' = r — R; to center 


dRu a Wn-n (R) = 

E! ZsZs-v^iRu-Ri-s'), (A16) 

Z — s' 

and 

9R lea R llsla ,W n - n (R) = 

- Z S Z S ,V%1 (R Is ~ R Va>) (1 - Su '<M 
+ E^, a „ z sZs"V^ a , (R ls - R i»s")5w5 S s'. (A17) 

This, leads to 



^s^s' r iq.Ri, 

x (R is - Ri's') Va (qA|s), 


(A18) 


and, finally, that 
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Z K e' 


(q-Rla+q'-Rl's') 


^( 2 q ,)(q-v)^-n (R)=EE -,E V/ W ^ (qA|s) (q ' A ' |S,) [- Z -' W i a 2' ( Ris - 

^ Isa l's'a' Wy/MsMs'UqxUq'x L 

,( 2 ) /"o. TT 


+ 5Z z „ a „ (R; s - R 1»,») 6u'6 sa . \ . (A19) 


r 


Appendix B: An iterative expression for arbitrary 
n—th order derivatives of V sc f with respect to the 
nuclear displacements 

One of the ingredients that are most difficult to cal¬ 
culate within an ab-initio framework is the higher or¬ 
der derivatives of V sc f (R, r). While the first and second 
order are needed to calculate the lowest order electron- 
phonon self-energies an extension of the theory to the 
regime of strong interaction requires the knowledge of an 
arbitrary order derivative 


In order to rewrite in terms of T^ n B we need to 
evaluate 


Pr ^|mp(R))r^-i 


mn 

Pk' 


■»-!) 

■Rn 


(Primp (R)>) T^:. : 

+ |mp(R)) (PiT^ 


P k' 

(ra-l) 


(B8) 


P k' 


p( n ) 


( r ) — 


d n V scf (R, r) 
<9Ri • • • <9R„ ‘ 


(Bl) 


In order to derive a close expression for T^ we use the 
following property of the nuclear momentum operator P/ 


Pi ,V scf (R, r) 




dV sc f (R, r) 


5Ri 


This identity can be iterated to give: 

P 2 , V sc f (R, r) 


Pi, 


More generally we have that 


d 2 V scf (R,r) 


<9Ri<9R 2 


T 


(n) 


RlR2"-Rn 


■p rp{n—l) 

^1’ 1 R2-Rn 


(B2) 


(B3) 


(B4) 


In order to evaluate the T^ matrix elements we start 
by noticing that, within DFPT, 




Now we define 

(n) 


T 


Rl ■ ■‘Rn 


^mp ^nk 


< = (nk|rE l) ..R„ ( r ) K k '), 


nn 

kk' 


(B5) 


(B6) 


and we plug into Eq. (1B41) a complete set of eigenstates 
obtaining 


Last step is, then, to use Eq. 


T 


(n) 


Ri"-R‘ 
(n-l) 


nn 

kk' 


T, 


Ro • "Rn 


'= E 

mp/nk 

(i) 


T, 


(i) 


to get 


rp{n—l) 

nm i?2' "Rn 

kp 


mn 

pk' 


T 

nm R 

kp 


P k " 


^n'k' 


^mp ^nk 


p* (Erl 


nn 

kk', 


(B9) 


In order to check the soundness of this approach we apply 
it to the first order case: 


t: 


(i) 


T 


(i) 


,= V ■ 

nn Z —j 

kk' mp/nk 


Ri 


T(°)| ; 

nm | mn 

kp pk' 


^mp ^nk 


T 0| 


E ■ 

mp^nk 


'kp 


r; 


(i) 


Ri 


mn 

pk' 


rp( 0) 


^n'k' 


nn 

kk' / 


But it is easy to verify that 

dV ac f (R, r) d \T e + V scf (R, r)] 


(9Ri 

which implies that 


0Ri 


(BIO) 


(BH) 


r (n) 

1 Ri-’-Rr, 


nn 

kk' 


* E 

mp^nk 


(nk (R) |Pj I |mp(R)>T<£-£ 


mn 

pk' 


- T, 


(n-l) 


R2 —Rn 


nm 

kp 


(mp|Pi|n'k') 


• (B7) 


ji(0) 


f — ^n'k'^n'ra^k'p • 


pk' 


(B12) 


Eq. (lB12l) .when used in Eq. (1B10I) . gives an identity 
that confirms the correctness of the iterative equation, 
Eq. (1B9I) . 

When Eq. (IB9I) is applied to the second-order deriva¬ 
tive it provides a well-known relation that connects the 
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second to the first derivativ e 21 ! 45 . 


rr i(2) 

1 RiR 2 


nn 

kk 


E 

rap 


T$ 

nm R 2 

_l_ T ^ 
ran 1 i ?2 

rpi}-) 
nm Ri 

ran 


k P 

P k 

kp 

pk 


^mp ^nk 


-4 lR2 ^k(R). (B13) 


More in general Eq. (IB9I) can efficiently rewrite the n- 
th order derivative in terms of the first order derivative 
and of the n-th order derivative of the electronic ener¬ 
gies, . But these derivatives can be efficiently 

calculated using DFPT or finite differences^^. 
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